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We present first results on the dynamics of monopoles in full QCD with Nf = 2 flavors of dynamical quarks. 
Among the quantities being studied are the monopole density and the monopole screening length, the static 
potential as well as the profile of the color electric flux tube. Furthermore, we derive the low-energy effective 
monopole action. 



1. INTRODUCTION 

In the dual superconductor picture of con- 
finement the crucial degrees of freedom are the 
color magnetic monopoles revealed after abelian 
projection. In the maximally abelian gauge [Q 
one finds that the string tension is almost en- 
tirely due to monopole currents 1^,^, and that 
the low-energy effective monopole action repro- 
duces the string tension and the low-lying glue- 
ball masses quite well Q. Furthermore, the 
monopole currents show non-trivial correlations 
with gauge invariant excitations of the vacuum, 
such as the action and topological charge den- 
sity The analysis of the monopole de- 
grees of freedom may thus provide important 
information about the confinement mechanism. 
Past studies have mainly been restricted to the 
quenched approximation. It will be interesting 
now to see how the vacuum reacts to the pres- 
ence of dynamical color electric charges. 

Our studies will be done on Nf = 2 dyna- 
mical gauge field configurations generated by 



the QGDSF and UKQGD coUaborations us- 
ing non-perturbatively 0{a) improved Wilson 
fermions [B|: 
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where Sp'' is the original Wilson action. 

The link variables U{s, ^) are brought into 
the maximally abelian gauge by maximizing the 
quantity ^ 

R = J2j2\^u{s,n)\' (1) 

with respect to gauge transformations g: 

Uis,fi)=g{s)U{s,fi)g{s + fi)-\ (2) 

using an SU{3) update of the simulated annealing 
algorithm first introduced in for SU{2). The 
abelian link variables are defined by 

with 

0i{s,fi) = arg[/ii(s,^) 
1 ^ 

- ■::^l^'rgUjj{s,^l)\nwd2TT- 
•^3 = 1 
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Figure 1. The monopole density for all clusters 
(x) and for large clusters (□). 



2. MONOPOLE DENSITY 

The first quantity we looked at is the monopole 
density. This is defined by 



(5) 



where ki{s,fi) are the monopole currents, which 
are obtained from the angles 9i{s,fi) and 
V is the lattice volume. The currents are sub- 
jected to the constraint '}Ziki{s, ^jl) = 0. The 
monopole currents fall into clusters of different 
lengths: large clusters are responsible for quark 
confinement, while short clusters are mainly due 
to ultraviolet fluctuations and have no effect on 
the conflning forces As large clusters we 

count the largest cluster of the system and all 
clusters that wind around the lattice at least once. 
It was shown Q (in pure SU{2) gauge theory) 
that the density of monopoles belonging to the 
large clusters is finite in the continuum limit. 

In Fig. 1 we show the monopole density for 
all clusters, and separately for the large clusters, 
as a function of (mjr/TOp)^. Also shown is the 
quenched result {mT^/mp = 1). The quenched 
data set, /3 = 6.0, was chosen to match the lattice 
spacing of the dynamical configurations, a « 0.1 
fm. The lattice volume varies between 16^32 
and 24'^48. We see that the monopole density 
is strongly affected by the presence of dynamical 
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Figure 2. The static potential (x) and the abelian 
static potential (□), with the self-energy sub- 
tracted. The upper figure is for [3 = 5.29, n = 
0.135, the lower figure for (3 — 6.0 in the quenched 
theory. 

quarks. On the dynamical configurations the den- 
sity is twice as large as on the quenched configu- 
rations with increasing tendency towards smaller 
quark masses. 

3. STATIC POTENTIAL 

The next quantity we have looked at is the 
static potential V{r). This is computed from 
smeared Wilson loops. We parameterize the po- 
tential by 



V{r) = - - + ^ 



(6) 



We also have looked at the abelian potential 
built from abelian link variables (p|). This de- 
fines an abelian string tension a . In Fig. 2 
we show V{r) and its abelian counterpart, and 
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Figure 3. The abelian static potential (x) divided 
into monopole (□) and photon (O) contributions. 
The upper figure is for (3 = 5.29, k — 0.135, the 
lower figure for (3 — 6.0 in the quenched theory. 

compare it with the quenched resuh. A fit gives 
0^"° ja — 0.94(8) for the dynamical configuration 
shown and a'^^ /a = 1.03(15) for the quenched 
case. 

It is instructive to separate the abehan gauge 
potential into a monopole and 'photon' part, 
6»(s, = 6''"™(s, /i) + 6'P''(s, m) where 



s')V~m^,{s'). (7) 



Here D{s) is the Coulomb propagator and mf^,j{s) 
is the number of Dirac strings piercing the abelian 
plaquette u(s, n, v). In Fig. 3 we show the abelian 
potential constructed from the monopole and 
'photon' part of the abelian gauge potential, re- 
spectively. We see that most of the string tension 
is due to the monopole currents, while the 'pho- 
ton' part of the potential is Coulomb-like, exactly 



1.35 ^ 
1.3 

IMO 

^ 1.25 r 
^ 1.2 |i 
1.15 
1.1 



10 

N 



15 



20 



Figure 4. The abelian string tension versus the 
number of gauge copies (x), and for iterative 
gauge fixing (□). 

as in the quenched theory. 

The quantity R may have many local maxima. 
Numerically it is difficult to get to the absolute 
maximum, which constitutes a problem, the Gri- 
bov problem. In this study we considered a sin- 
gle gauge copy per configuration. To estimate 
the error arising from incomplete gauge fixing we 
created a sequence of random gauge copies on 
each of our configurations at /3 = 5.29, k = 0.135, 
following 1^. In Fig. 4 we show the abelian string 
tension as a function of the number of gauge 
copies considered (thereby always computing a^^ 
on the copy with the largest value of R). We also 
show the result for iterative gauge fixing. We see 
that the bias induced by incomplete gauge fix- 
ing is about 6%, while iterative gauge fixing may 
give results that are wrong by 0(20%). Therefore 
careful gauge fixing is necessary. Similar results 
hold for other observables. 

4. SCREENING LENGTH 

The screening length ^ of a plasma of 
monopoles is defined by the exponential decrease 
of the magnetic flux through a sphere of radius r 
around the monopole: 

$(r) = $oe"''^^. (8) 

A screened plasma of static monopoles is known 
to generate a string tension ||] 

acxpe (9) 
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If this is true for dynamical monopoles as well, 
we would expect the screening length to be signifi- 
cantly smaller on dynamical configurations. Con- 
sidering the monopoles in the large clusters only, 
we find 
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We see that the screening length is indeed 0(40%) 
lower in the dynamical case. We interprete this, 
as well as the increase in the monopole density, 
as pairing effect jlj] . 

5. ABELIAN FLUX TUBE 

The chromoclectric flux tube between static 
color electric charges forms an interface between 
the superconducting (confining) phase and the 
conducting phase, and in this capacity bears valu- 
able information about the properties of the QCD 
vacuum. Because we are primarily interested 
in the dynamics of monopoles, we may restrict 
ourselves to abelian flux tubes being excited by 
abelian Wilson loops. The abelian Wilson loop is 
defined by 



W^(i?, T) = i Tr n m(s, m) = ^ Tr e''^ , 



(10) 



sec 



where C is a rectangular loop of extension RxT, 
9c = diag[6'i, 6*2, 6'3]c. In this section we shall 
study the spatial profile of the flux tube. This 
is done by looking at correlations of appropriate 
operators O = diag[C'i, C's] with the abelian 
Wilson loop operator: 



l (Tr0(g)Tre''^^) 1 

3 (Tre'^c) 3^ ^ 

for C-parity even operators, and 

(Tr (C)(s) e'^^O) 



(11) 
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Figure 5. The action density PA{s)rQ of the 
abelian flux tube in full QCD (mjr/mp « 0.6) and 
quenched QCD. 



for C-parity odd operators, where s is taken to lie 
on the i? X T plane at T/2. 

The first quantities we looked at are the action 
density 



TrO(s) = I XI cose^{s,p,i^). 
and the color electric field 



0,(s)=i?j(s) = i^,(s,4,j), 
where 

0,{s,p,iy) = -i\n{e'''^''^^''''>). 



(13) 



(14) 



(15) 



In Figs. 5 and 6 we show our results on a 16^32 
lattice, where we took R = 10, which corresponds 
to a separation of « 1 fm, and T = 6. The action 
density is somewhat larger in the dynamical case, 
while we see little difference in the distribution of 
the electric field between full and quenched QCD. 
A fit of Eu — const. exp(— r^/A), restricting r± 
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Figure 6. Distribution of the color electric field. 

to r± > 0.25 fm, gives the penetration length 
A = 0.15(2) fm in fuU QCD and 0.17(1) fm in the 
quenched theory. 

Next we looked at the local monopolc density 

TrO(s) = ^^|fc,(s,/i)|. (16) 

The result is shown in Fig. 7. In this case 
(1/3) (TrO) was not subtracted in eq. (|lj). This 
figure shows again that the monopole density, 
outside the flux tube, is more than twice as 
high on the dynamical configurations as in the 
quenched case. Inside the flux tube the monopole 
currents are strongly suppressed, as expected. 
Due to the high monopole density outside the flux 
tube, the interface appears to be much 'harder' in 
fuU QCD. 

We expect the monopole current to form a 
solenoidal supercurrent which, similar to a coil 
in electrodynamics, constricts the color electric 
fields into flux tubes. This is expressed by the 
dual Ampere law 

k = V X E. (17) 




Fi gure 7. The local monopole density p{s^ Tq in 
the vicinity of the ahelian flux tube. 

We find that the monopole current is indeed az- 
imuthal, as can be seen in Fig. 8. In Fig. 9 we 
compare the l.h.s of eq. ( p7| ) with the r.h.s. and 
find that Ampere's law is very well satisfied on our 
dynamical configurations. Ampere's law was al- 
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Figure 8. The solenoidal monopole current in a 
plane perpendicular to the abelian flux-tube, posi- 
tioned at R/2. 
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Figure 9. Comparing k and V x E. 

ready verified in pure SU(2) gauge theory in jl5|. 

We expect the flux tube (string) to break even- 
tually if the static charges are pulled apart be- 
yond a certain distance. That distance is ex- 
pected to be between 1.5 and 2 fm, depending 
on the mass of the quark. In Fig. 10 we show 
the distribution of the color electric field in the 
vicinity of the flux tube on the 24'^48 lattice 
at f3 — 5.29, K — 0.1355, corresponding to an 
TTiTr/rnp ratio of « 0.7, for R — 22 which amounts 
to a separation of « 2 fm. We see no sign of 
string breaking yet. We notice though that the 
electric field vectors show some level of noise half 
way between the charges. 

The penetration length and half width of the 
flux tube turn out to be « 0.16 fm and « 0.3 
fm, respectively. It is striking that the abelian 
flux tube does not show any broadening towards 
larger distances from the static charges, as one 




Figure 10. Distribution of the color electric field 
on the 24^48 lattice for R ^ 22,T = 5 in full 
QCD. 

would expect if the fluctuations of the string 
were described effectively by the Nambu-Goto ac- 
tion |l^. (This effect has also been observed in 
pure SU{2) gauge theory jl^.) This could mean 
that the abelian string is described by the Ra- 
mond string action instead, which, for example, 
does not show such a broadening effect [n8| . 



6. EFFECTIVE MONOPOLE ACTION 

We have seen that the vacuum undergoes sev- 
eral changes if dynamical color electric charges 
are introduced. We shall study now how this will 
affect the effective monopole action. 

There are three types of monopole currents, of 
which two are independent. For simplicity we 
take into account only one of them, thus integrat- 
ing out the other two For the time being, we 
assume the form of the effective monopole action 
in full QCD to be the same as in the quenched 
theory ||2^. This is composed of 27 types of two- 
point interactions, one four-point interaction and 
one six-point interaction: 



29 



(18) 



where d are the coupling constants which need 
to be determined. This we will do by employing 
an extended Swendsen method Q. Explicitly, 
we have: 



7 



Two-point interaction {i — 1, . . . , 27) 

4 

s /i— 1 

4 

S'2(fc) = ^^fc(s,^)fc(s + /i,/i), 

S fl—1 

Ssik) = ^^fc(s,Ai)fc(s + :>,Ai); (19) 



Four-point interaction 



s \ai=-4 / 



(20) 



Six-point interaction 

529(fc) = Ef E ^(^'^)') ■ (21) 

s \ai=-4 / 

Our calculations are done on 0(50) gauge fixed 
configurations each on the 24^48 lattice at /3 = 
5.29, K = 0.1355 and on the 16^32 lattice at /3 = 
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Figure 11. The coupling constant of the monopole 
self-interaction, Gi, as a function of the physical 
length scale b. The symbols are: n — 1 (9) , n — 2 
(U), n = 3 (♦), n = 4 (A) for full QCD, and 
n = 1 (O), n = 2 (□), n = 3 (0), n = 4 (A) for 
quenched QCD. 
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Figure 12. The coupling constants G2 and G3 as 
a function ofb. The symbols are as in Fig. 11. 

5.29, K = 0.1350 [1). For comparison, we have 
also done calculations in the quenched theory on 
24^48 lattices at /3 = 5.8, 6.0 and 6.2, and on a 
16^32 lattice at /3 = 6.0. 

We have employed a type-II block spin trans- 
formation [P2| with up to n — 1 = 3 blocking 



steps. In Fig. nil we show the sell-coupling Gi as 
a lunction of the physical length scale b — na for 
both full and quenched QCD. We see that Gf^'' 
is systematically smaller than q'^^^^^^'^ fQp &, 
in agreement with the higher monopole density 
in full QCD. A necessary condition for monopole 
condensation is Gi < In 7. This is achieved for 
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Figure 13. The coupling constants G28 and G29 of 
the four-point and six-point interactions, respec- 
tively, as a function of b. The symbols are as in 
Fig. 11. 

b > 0.27 fm in full QCD and for b > 0.35 fm in 
the quenched theory. In Fig. |2| we show the cou- 
pling constants G2 and G3. We see that Gj"'' < 



G 



quenched 



as in the previous case, while G3 shows 



the opposite behavior: G|"" > G(3)^"'="'='^'='^. For 
the other coupling constants we find little differ- 
ence between full QCD and the quenched theory. 
In Fig. |l^ we show, as an example, G28 and G29. 
Note that the self-coupling is the dominant cou- 
pling. 

To shed some more light on the dynamics of 
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Figure 14. The two-point coupling of monopole 
currents fc(s,/i) and k{s',fi) as a function of dis- 
tance d = \l^n{sii ~ s'Y /i*^^ and quenched 



QCD. The symbols are: full QCD quenched 
QCD (■) and Coulomb propagator (0). 

the monopoles we have looked at the coupling G 
of two (n = 1) parallel monopole currents, fc(s, /i) 
and fc(s',/x), as a function of their distance d = 

y^^^(sp — s'^y. In Fig. |lj we show G together 

with the Coulomb propagator. We see that at 
distances d ^ 2 the interaction becomes weaker 
than Coulomb in both full and quenched QCD. 
This is a result of the screening effect discussed 
earlier on. It comes as a surprise though that 
we see no difference between full QCD and the 



9 



quenched theory. We would have expected the 
screening effect to be considerably stronger in full 
QCD. 

Our task for the future is to check for scaling of 
the effective monopole action, as this has already 
successfully been done in the quenched theory. 

Furthermore, it would be useful to cast the ac- 
tion into the form 

S{k) — 5coulomb + "Ssclf + '5'4 — point + <5'6-point,(22) 

to make contact with the dual Ginzburg-Landau 
theory, and the string model, and determine their 
parameters. This would allow us to do truly 
quantitative, analytic calculations p^]. 

7. CONCLUSIONS 

We have had a first look at the dynamics of 
color magnetic monopoles in full QCD. We found 
striking differences between full QCD and the 
quenched theory. 

How can this be understood? Monopoles are 
at least partly induced by instantons The 
fermion determinant causes instantons and anti- 
instantons to attract each other, with a force that 
increases with decreasing quark mass. The effect 
is that isolated instantons are suppressed, and in- 
stantons and anti-instantons are forced to form 
(overlapping) pairs. This will increase the density 
of (anti-)instantons, and consequently the density 
of monopoles. 
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